Extended dynamic Mott-transition in the two-band Hubbard model out of equilibrium 
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We reformulate the time-dependent Gutzwiller approximation by Schiro and Fabrizio [PRL 105, 
076401 (2010)] in the framework of slave-boson mean-field theory, which is used to investigate the 
dynamical Mott transition of the generic two-band Hubbard model at half filling upon an interaction 
quench. Inter-orbital fluctuations lead to notable changes with respect to the single-band case. 
The singular dynamical transition is replaced by a broad regime of long-lived fluctuations between 
metallic and insulating states, accompanied by intriguing precursor behavior. A mapping to a spin 
model proves helpful to analyze the different regions in terms of the evolution of an Ising-like order 
parameter. Contrary to the static case, singlet occupations remain vital in the Mott-insulating 
regime with finite Hund's exchange. 
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Introduction. — Enormous advances in the physics of 
ultracold gases [l|, |2j have evoked strong interest in time- 
dependent (TD) phenomena, which can be provoked and 
studied in trapped cold atom systems without all compli- 
cations that instead arise in solid-state materials. This 
gives the unique opportunity to investigate physical real- 
izations of prototypical models for interacting particles, 
like Bose or Fermi Hubbard models, and examine fun- 
damental questions not only at equilibrium [3|, but also 
in out-of-equilibrium conditions [4-6]. The simplest pro- 
tocol to drive a system out of equilibrium is a sudden 
change of its Hamiltonian parameters, for instance an in- 
teraction quench. An intriguing issue thereof is the possi- 
ble trapping within metastable configurations that have 
no stable equilibrium counterpart |7H11|. Eventually, 
these metastable states decay at long times into some 
thermal configuration, i.e. thermalization occurs. The 
detailed influence of thermalization onto TD phenomena 
is a subject to its own [a. 1 12 - 1 15] . Considering the simple 
Fermi Hubbard model, most investigations concentrated 
so far on the single-orbital case |l2i . Il6] . However materi- 
als in nature are most often ruled by multiorbital degrees 
of freedom. To what extent multiorbital processes mod- 
ify the encountered single-orbital Hubbard physics out of 
equilibrium is a highly-relevant question. 

In this work we report about qualitative new physics 
close to the dynamic Mott transition of the canonical 
two-band Hubbard model upon an interaction quench. 
The singular dynamic Mott transition in the single-band 
case is smeared out to a broad region with long-time fluc- 
tuations between metal and insulator. Unique precursor 
behavior takes place at the borders of this region. 

Theoretical approach and model. — The dynamic Hub- 
bard model is solved within TD slave-boson mean- 
field theory (SBMFT) that merges the rotational in- 
variant equilibrium slave-boson method [171 ] with a non- 
equilibrium scheme recently proposed by Schiro and Fab- 
rizio [lg, [l8|, [l9( via the Gutzwiller representation [20( • 



In the latter scope the time evolution is described via 
coupled first-order differential equations for Slater deter- 
minants ipo and Gutzwiller projectors (f>, i.e. 
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whereby H [(/}] denotes the renormalized free Hamiltonian 
and TL loc the local interacting one (note that explicit site 
dependance is omitted). To reexpress these equations in 
SBMFT an unitary transformation of the Slater deter- 
minants \ipo) into eigenstates i/ k of H[4>] in momentum 
(k) space is performed. Hence the former eigenstates 
are identified with the quasiparticle (QP) degrees of free- 
dom. For the complete transformation, utilizing a, b for 
the eigenvalue labeling and A, B, C for the local basis 
states, one obtains 
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with a, B denoting the respective orbital-spin combi- 
nation. The matrices R aa ' =R aa ' [4>] renormalize the 
free dispersion e^a and form the QP-weight matrix 
Z=RR t [l7J|. Equations |0 0} show the coupled time 
evolution of QP states z/ k and local slave-boson ampli- 
tudes 4>ab ■ Within mean-field versions of Gutzwiller and 
slave-boson techniques, those quantities are subject to 
certain constraints, reading (for each point in time t) 
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A numerical solution of eqns. (|3I4[) is achieved by using 
an adaptive Runge-Kutta scheme of order 5/6 [21] . 
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FIG. 1. (Color online) QP weight Z(t) and occupation prob- 
ability #p of local sectors with particles p— 0,1,2 for the inter- 
action quench to Uf =0.2. Dotted lines denote period of local 
states oscillations. The inset shows the Fourier transform 
(absolute value of Fourier coefficients in logarithmic units) of 
Z(t). See the text for use of the circles. 



We study the canonical two-band Hubbard Hamilto- 
nian H = H + J2i W° c with a nearest-neighbor hopping 
t that defines the kinetic part H . In detail it reads 
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where i, j are site indices, m, ml run over orbitals 1,2 and 
<T=t, 4- marks the spin projection. A three-dimensional 
simple-cubic dispersion is used and thus the parametriza- 
tion U'=U-2J, U"=U-3J proves adequate (H, [H|. 
The value of the hopping t is such that the half- 
bandwidth D is the energy unit. To investigate an in- 
teraction quench the initial t/j is set to zero and the 
local interaction varies in time as U(t)=Uf ®(t) and 
J{t)=qUf &(t) , i.e. jumps from zero to Uf and qUf 
at t=Q. In the following we focus on the paramagnetic 
half-filled scenario. 

Results. — First the limiting case J=0 (q=Q), i.e. 
U=U'=U", is examined. The quench leads to non- 
vanishing oscillations in the physical observables due to 
the lack of quantum fluctuations in the present formal- 
ism 0. In Fig.[T]the TD QP weight Z ma (t)=5 12 6 n Z (t) 
and occupation probability #p of local particle sectors is 
displayed for small Uf. While Z (t) and the occupation 
in the p=0,l sectors oscillate in phase, the two-particle 
sector commutes with a 7r phase shift thereto. As ex- 



pected, minima in Z (t) amount to maxima in =ff=2 and 
vice versa. The Fourier transform (inset in Fig.[T]) reveals 
the occurrence of two main frequencies w\, U2 as well as 
sidebands. This is in contrast to the single-band and the 
quasi-decoupled (U'=U"=Q) two-band case, where only 
a single frequency on the order of uj\ shows up. While 
the latter frequency in the fully-interacting case is sen- 
sitive to changes in the hopping t, U2 remains rather 
unaffected. Hence LU2 originates from interorbital pro- 
cesses and is absent for the quasi-decoupled-band case. 
The encircled frequencies are integer multiples of wi, W2 
or sidebands and appear due to Fourier transformation 
on a finite-time interval (note logarithmic units). The 
time interval was set to 7500 D^ 1 to provide suitable 
high-frequency resolution <(1000-D) -1 . 

The overall relevant-frequency behavior from the 
Fourier transform of Z(t) with increasing Uf is depicted 
in Fig. [2] with comparison to the U'=U"=0 case. For 
the latter, the bands are independent and we recover the 
results of the single-band study [16(. A dynamic Mott 
transition occurs at Uc =1.325 eV indicated by a loga- 
rithmic divergence in the single frequency wis, smeared 
out due to the numerical integration. With interorbital 
terms, multiple relevant frequencies appear with growing 
Uf. In fact, 1M2 consists of two coupled interorbital con- 
tributions, namely u>2top and tU2bot- Both contributions 
tend to zero close to the dynamic Mott criticality, which 




FIG. 2. (Color online) Top: main-frequency evolution with 
U f for U'=U"=0 and U'=U"=U from Fourier transforming 
Z(t). The dotted lines separate three different regions (see 
text and Fig. |3J. Bottom: Time-averaged Z as a function 
of Uf. Dashed lines mark the equilibrium behavior with the 
black dot denoting the associated Mott-critical point. Grey 
hashed area indicates precursor behavior. 
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FIG. 3. Level diagram for the active frequencies in the case 
U'=U"=U within the different regimes dictated by Uf. 



we believe to occur near [/ c (see below). Three obser- 
vations from Fig. [2] are vital. First the splitting of lo-2 
for 0.3<[//<1.0 which can be visualized within a level 
scheme (see Fig. [3]). That scheme shows the onset of 
splitting between levels connected by a>i when approach- 
ing the Mott critial regime. Second a broad region of 
mostly noisy Fourier spectra appears for 1.15<f/<2.375 
without a clear structuring. Albeit already insulating, a 
definite level structure then reemerges for [// >2.375. 

Whereas U2 saturates for large interaction-quenches, 
lo\ linearly rises with the same slope as in the quasi- 
decoupled-bands case. Figure [2] also shows the quantity 
Z=Z(T)=^ L dt Z, based on time-averaging over a long 
time interval T (not to be confused with a basic period 
of oscillation), in comparison with the equilibrium QP 
weight. If not specified T was set to 7500 D~ l . A well- 
defined Z is only accessible outside the broad [//-range 
with noisy Fourier spectrum. Interestingly, just before 
entering the latter regime a small range appears where Z 
of the U'=U"=0 model is smaller than in the one with 
including interorbital terms. Thus interorbital interac- 
tions allow to reduce the standard correlation measure 
in certain out-of-equilibrium cases. 

Better understanding of the dynamical transition at 
J=0 can be gained by an Ising-spin representation of 
the model that is fully equivalent to slave-bosons, and 
reads El 
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with spin S=2 and where r=6 is the lattice coordination 
number while — J=— 2/3 is the energy per site of the 
non-interacting ground state. Within mean-field, i.e. as- 
suming a variational wave function Yii \&{Siz)}, the spin 
model above becomes identical to the slave boson mean- 
field theory at half-filling if \${S iz )) = £p=o P \S tz = 
p — 2) , where <f> p is the original slave boson at site in par- 
ticle sector p. Because of half-filling, 4>, 



i- p , so that 
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FIG. 4. (Color online) Order parameter (S x ) for U'=U"=0 
(dashed) and U'=U"=U (full). Inset shows (S x ) (T) in the 
three different regions. The quasi-decoupled two-band case 
stabilizes only in regions I and III (divised by the vertical 
dashed line). Region II appears with interorbital terms. The 
precursor regime is visualized by the grey hashed area. 



(Siz)=0- Metallic coherence is signaled in the spin model 



by a finite Ising order parameter (Si X )=$t 
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while the inco- 
herent Mott insulator has (Si X )=(Sf z )=Q. The initial 
non-interacting state is characterized by (Si X )=2 and 
(Sf z )=l, i.e. energy per site E=— 2/3+C//2, conserved 
during the unitary evolution. If, like in the single-band 
case [16|, we assume the dynamical Mott transition to 
occur when the energy equates that of the Mott insula- 
tor, i.e. ^=0, then we would expect Ut\ =4/3, which 
is also equal to the single-band value for the dynamical 
Mott transition. In fact, this argument predicts one and 
the same value C/ c =4/3 for any A^-band simple, i.e. J=0, 
Hubbard model at half-filling, indeed close to the numer- 
ical value. 

Figure[4]exhibits the time-averaged quantity (S x ) with 
increasing Uf, indeed classifying three different regions. 
Hence (S x ) does serve here as an order parameter, where 
Regime I is metallic with a finite (S x ). However there is 
a crossover of the contributions Mi, M 2 for U'=U"=U, 
showing that correlations between local states of the one- 
(three-) and two-particle sectors are less affected by in- 
creasing Uf than those between zero-(four-) and one- 
(three-) particle sectors. For U'=U"=Q the order param- 
eter jumps to zero at Uc , signaling the dynamic metal- 
insulator transition into a regime III. With finite interor- 
bital interactions the behavior is surprisingly more in- 
triguing. At Uf =1.15 a precursor regime starts with in- 
creased fluctuations of growing frequency. The integrated 
components Mi, M 2 are specifically sensitive thereto. 
This may be connected to the breakdown of U2 shown 
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FIG. 5. (Color online) Time-averaged diagonal renormaliza- 
tion R and QP weight Z for U/J=5. Black dot marks the 
Mott-critical point of the associated equilibrium model. In- 
sets show R(T) in the three different regions. 




FIG. 6. (Color online) Time-averaged multiplet occupation 
probabilities for U/J=5 with the equilibrium probabilities 
(dashed lines) for comparison. The Mott-critical point of 
the associated equilibrium model is described by a black dot. 
Insets reveal explicit time evolution of multiplet occupation 
probabilities in different regions. 



in Fig. [21 leading to a noncontinuous evolution of corre- 
lations between the local particle sectors. In the region 
1.375<Uf<2.0 the order parameter then indeed vanishes, 
but only after an extremely long time. This extended 
Mott-insulating transition (EDMT) is characterized by 
a chaotic-like time evolution of the QP weight and the 
multiplet-occupation probabilities. This is in accordance 
with the noisy Fourier spectra. Note that the EDMT 
is a unique feature of the general two-band model with 
interorbital interactions. The static mott-critical interac- 
tion strongly differs if the orbitals are explicitly coupled 
(U c ^2.65) or not (Z7 C ~4.0). In contrast the boundary 
of the dynamic metallic region I is quite independant of 
the coupling type, but a qualitative difference is intro- 
duced via the appearance of an intermediate region for 
interaction-coupled orbitals. 

The EDMT is not an artifact of the ( J=0) case but 
a general feature of the general dynamic two-band Hub- 
bard model. Since (S x ) is proportional to the R ma- 
trix (see eq. Q), the quantity R serves as a suitable 
order parameter in the original model and is depicted for 
U/J=l/q=5 (including all interaction terms in eq. (6)) in 
Fig. The precursor regimes as well as the occurrence of 
the EDMT persists. In fact with finite J this intermedi- 
ate region is even broadened (l.I25<[//<2.225) in the R- 
Uf diagram compared to the J=0 case (1.375<£//<2.0). 
In addition, Fig. [5] shows the evolution of the time- 
averaged multiplet-occupation probabilities for finite J. 
For small Uf they are close to the values in the equilib- 
rium model, but lack the strong polarization in the in- 
sulating region. While in the static case only the triplet 
channels survive the Mott transition, the dynamic Mott 
region III displays also finite singlet fillings. 



Summary. — We have investigated an interaction 
quench in the canonical two-band Hubbard model with 
the TD-SBMFT scheme allowing for complete rotational- 
invariance. The limited case with sole intraorbital inter- 
action terms leads to a quasi-decoupled two-band model 
with dynamic characteristics reminiscent of prior single- 
band studies |16| . On the other hand, when introducing 
the relevant interorbital interactions, novel physics ap- 
pears out of equilibrium. An intermediate region with 
long-time chaotic-like fluctuations in the physical ampli- 
tudes and high-frequency metal-to-insulator fluctuations 
emerges and replaces the singular Mott-transition point. 
This replacement goes along with precursor regimes ex- 
hibiting high-frequency fluctuations in time. Finally the 
appearance of chaotic-like behavior may be connected to 
the non-integrable classical coupled-pendulum problem, 
which indeed displays chaotic orbits. Extensions beyond 
mean-field are needed to reveal, if such analogies hold in 
the complete quantum-fluctuating scenario. 
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